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Abstract 

We present high-resolution (0."2) near- infrared images of the area surrounding the 
Herbig Be star LkHa 234 taken with the Coronagraphic Imager with Adaptive Optics 
(CIAO) and the adaptive optics on the Subaru Telescope. The near-infrared (J, H, K, 
V and M' bands) images reveal circumstellar structures around LkHa 234 in detail. 
Eight young stellar object (YSO) candidates (Object B, C, D, E, F, G, NW1, and 
NW2) were detected at 2"-ll" from LkHa 234. Objects B and C are likely variable 
stars, which is consistent with a young evolutionary status. Three objects (LkHa 234, 
NW1, and NW2) were identified in the 11.8 /im (SiC) and 17.65 fim images obtained 
with the Long Wavelength Spectrometer (LWS) on the Keck Telescope. NW1 and 
NW2 are thought to be embedded young stars. We suggest that NW1 and not 
LkHa 234 is the source illuminating the reflection nebula west of LkHa 234, although 
Object G may be another candidate. In our images, these objects are located at 
the center of the 2 /im polarization, and NW1 resides at the center of a monopolar 
cavity. The SED of NW1 suggests that it is a YSO with a spectral type of B6-7. 
Object F and G were discovered by our observations at l'/9 and 2'/3, respectively, 
from LkHa 234, and their proximity to LkHa 234 suggests that they may be its 
lower-mass companions. 

Key words: infrared; stars — ISM; reflection nebula — stars; binaries; general — 
stars; pre-main sequence 



Based in part on data collected at Subaru Telescope, which is operated by the National Astronomical 
Observatory of Japan. 
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1. Introduction 



To understand star and planet formation mechanisms, it is necessary to understand 
the properties of circumstellar environments such as cluster density, circumstellar disks and 
envelopes. Although many stars are formed within clusters (Lada & Lada 2003), the process of 
cluster formation is not clearly understood. Other unanswered questions remain; for example, 
what is the difference between the formation of single and binary components? Are companion 
stars formed via fragmentation of a molecular core or disks (Wheelwright et al. 2010, Baines et 
al. 2006)? Moreover, the nature of protoplanetary disks affects the process of planet formation. 
Past observations at radio and infrared (IR) wavelengths show that circumstellar disks are 
ubiquitous around young stars (e.g., Williams & Andrew 2006). Observations at high- angular 
resolutions have frequently revealed a complex distribution of circumstellar material, which 
might be attributed to stellar multiplicity and properties of the local star-forming environments 
such as initial disk masses (e.g., Fukagawa et al. 2010). Therefore, it is important to investigate 
the environments of young stellar objects (YSOs) in the phase of star formation on a small 
spatial scale. 

Herbig Ae/Be (2 M < M* < 10 M & ) stars are intermediate-mass pre-main-sequence 
stars that are in evolutionary stages similar to lower mass T Tauri (0.1 M & < M* < 2 M©) 
stars. It is thought that the formation mechanism of low-mass stars is by the gravitational 
collapse of a single cloud core (Palla & Stahler 1993). In contrast, the formation and evolution 
mechanisms of massive stars (M* > 10 M©) are still under debate (Palla & Stahler 1993, Wang 
& Looney 2007, Zinnecker & Yorke 2007). Because they reside, at the boundary between low- 
and high-mass stars, Herbig Ae/Be stars can provide significant insight into star formation. In 
addition, they can be observed more easily than massive stars because their environments are 
comparatively simple (Testi et al. 1999, Wang & Looney 2007, Lada & Lada 2003, Hillenbrand 
1997). Further, because Herbig Ae/Be stars are more common than massive stars, a large num- 
ber of stars including much closer ones can be investigated. However, the detailed circumstellar 
environments around Herbig Ae/Be stars, especially Herbig Be stars, remain less understood 
than those around T Tauri stars. This is partly because it is difficult to detect such stars owing 
to their rapid evolution. Moreover, because the youngest Herbig Ae/Be stars are typically 
located at large distances (~kpc), high-spatial resolution is key to understanding their nature. 

The Herbig Be star, LkHa 234, is located in the NGC 7129 star-forming region at 
1.25 kpc (Shevchenko & Yakubov 1989). LkHa 234 is a star of spectral type B5~7, mass 8.5 M©, 
luminosity ~1700 L©, and estimated age ~10 5 yr (Hillenbrand et al. 1992, Fuente et al. 2001). 
The NGC 7129 star-forming region was previously observed at various wavelengths from the 
optical to radio. Another deeply embedded YSO located at 3" from LkHa 234 was detected by 
2 /xm polarization observations; this could be a protostar that illuminates the reflection nebula 
(Weintraub et al. 1994, Weintraub et al. 1996). Millimeter- wave interferometric observations 
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reveal two objects (IRS 6 and MM1) at a few arcseconds from LkHa 234, and a large CO 
outflow is observed near them (Fuente et al. 2001). The position of the protostar inferred by 
the 2 /im polarimetry is different from that inferred from the millimeter-wave interferometric 
observations. Other clumps (at 4"4, 3"0, 2''1 and 2''2 from LkHa 234) and H2O maser spots 
were detected by observations of the radio continuum and the H2O maser line with the Very 
Large Array (VLA; Trinidad et al. 2004, Umemoto et al. 2002). These previous observations 
showed that the CO outflow detected around this region is driven not by LkHa 234 but by other 
objects. In addition, H2 and [S II] jets observed in this region do not point toward LkHa 234, 
which suggests that there are other unseen sources producing the jets (Ray et al. 1990, Cabrit 
et al. 1997). To understand the unresolved sources driving these complicated outflows, it is 
essential to investigate the masses and locations of objects around LkHa 234. 

In this paper, we report observations of the Herbig Be star LkHa 234 using the stellar 
coronagraphic camera CIAO on the Subaru Telescope and the LWS on the Keck Telescope. 
Our high-resolution images reveal detailed circumstellar structures around LkHa 234 and eight 
YSOs at 2"-ll" from LkHa 234. 

2. Observations 

2.1. Near-infrared Observations 

LkH« 234 was observed using the adaptive optics (AO) with Subaru/CIAO (Tamura et 
al. 2000, Murakawa et al. 2004). Near-infrared images were obtained in the J (1.25 /im), H 
(1.64 /im), and K (2.20 /im) bands on July 15, 2006. In addition, the images were taken in 
the J band alone on May 25, 2008, and the V (3.77 /im) and M' (4.68 /im) bands on July 
10, 2008. The sky was clear all three nights. CIAO utilizes a 1024 x 1024 InSb ALADDIN 
II array. The pixel scale of the medium- resolution mode was 0''0213 pixel -1 , and the field of 
view was 21: "8 x 21."8. The observations (bands, exposure times, number of coadds, size of the 
occulting masks, photometric standard stars, and natural seeing) are summarized in Table 1. 

Adaptive optics (AO) allowed us to achieve high-spatial resolution given by the full- 
width at half-maximum (FWHM). On July 15, 2006, a spatial resolution was 0"15-0"2 in the 
J, H, and K bands, and on May 25, 2008, a spatial resolution of 0"17 was achieved in the J 
band. Finally, on July 10, 2008, spatial resolutions of 0''22 and 0''20 were achieved in the L' 
and M' bands, respectively. 

2.2. Mid-infrared Observations 

Mid-infrared observations were obtained with the LWS (Jones & Puetter 1993) on the 
W. M. Keck I Telescope. We obtained images using the 11.8 /im "SiC" filter (AA = 2.32 /mi) 
on July 26, 2004 and using the 17.65 /im filter ( AA = 0.85 /im) on August 28, 2004. Conditions 
were clear and dry on both nights. As measured from the Caltech Submillimeter Observatory 
(CSO) radiometer, the precipitable water vapor was 2.2 mm during the 11.8 /im SiC observa- 
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tions and only 1 mm during the 17.65 /im observations, thereby providing good atmospheric 
transmissivity. The LWS was used in the chop-nod imaging mode with a 10" throw at a position 
angle of degrees and a chop frequency of 2.5 Hz. The total on-source integration time was 
285 s for observations with each filter. The point-spread function (PSF) FWHMs were 0."32 at 
11.8 /im and 0"47 at 17.65 /im. For photometric calibration, we observed HD 197989, a mid- 
infrared photometric standard, from the catalog of Cohen et al. (1999) in a similar manner. 
The observations (bands, exposure times, photometric standard stars, and natural seeing) are 
summarized in Table 1. 

3. Data Reduction 

3.1. Near-infrared Data Reduction 

The images were processed using the Image Reduction and Analysis Facility (IRAF) 
packages (Tody 1986) 1 for dark subtraction, flat-fielding with sky flats, sky subtraction, and 
bad pixel correction. For the V and M'-band images, sky emission was cancelled by subtracting 
each pair obtained, by nodding the telescope in the north-south direction by 10". The total 
exposure time was just half of the total time given above for Objects B, D and E (Section 4 for 
designations of objects), each of which were covered by only one pointing of the nodding pair. 

Nine objects were detected in our near- infrared images in the field of view of 21. "8 x 21. "8 
(LkHa 234, Objects B-G, and NW1-2; Figures 1 and 2). We conducted PSF fitting and aper- 
ture photometry using the DAOPHOT and APPHOT IRAF packages. PSF fitting photometry 
was applied for Objects F and G in the JHK bands with an empirical PSF template made 
from Objects B and C, because the magnitudes for these objects could not be accurately mea- 
sured owing to contamination from the skirts of the PSF of LkHa 234. Object B, C, D, and 
F were measured by aperture photometry. The photometric values were corrected using the 
measurements of a photometric standard star with the same aperture (of about four times 
the FWHM). Relative photometry to Object B was also performed with the same aperture 
as Object B, which produced results consistent with the PSF fitting and aperture photometry 
results. Photometry in the V and M' bands was performed by aperture photometry. 

3.2. Mid-Infrared Data Reduction 

The LWS data were reduced using a custom Interactive Data Language (IDL) pipeline 
based on the LWSCOADD routines made available by Keck Observatory plus extensions for 
atmospheric transmission correction, photometric calibration, and source detection. First, each 
chop-nod set was reduced by standard double-differencing to remove background emission. 
No flat fielding was performed. Bad pixels were identified and replaced by the median of 

1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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neighboring pixels. The data were then rotated using a cubic spline interpolation to orient 
north up. All individual chop-nod-subtracted images were then registered via Fourier cross- 
correlation and summed into a mosaic image. 

To estimate a correction for atmospheric transmission losses, we first calculated a grid of 
Mauna Kea atmosphere models within a range of airmasses and water vapor contents, using the 
BFATS atmospheric transmission code (Roe 2002). For each model, we computed an overall 
atmospheric transmission weighted by the relevant transmission profile for each filter. Then, 
for each observation, given the measured water vapor (from the CSO radiometer) and airmass, 
we obtained a transmission correction by linearly interpolating between the closest models in 
the precomputed atmosphere grid. 

Observations of standard stars were reduced in the same way as the science data. We 
followed the method provided by Cohen et al. (1999) for flux calibration based on determination 
of the instrumental relative spectral response (RSR) consisting of the product of transmission 
curves for the atmosphere, instrumental optics, filter, and detector quantum efficiency We 
also derived aperture correction factors based on the radial profile of the standard star's PSF. 
Finally, we evaluated the absolute precision of these flux calibrations by comparing with the 
results of other authors; LWS photometry was calibrated using a near-identical procedure by 
McCabe et al. (2006), and the agreement in absolute flux calibrations is within a few percent. 

4. Results &z, Discussions 

4.1. Nebulae Around LkHa 234 

Our imaging observations reveal remarkable nebulae around LkHa 234. One of them, 
located at approximately 3" (~4000 AU) west of LkHa 234, extends more than 10" westward. 
In addition to appearing in the J, H, and K bands, it appears in the V and M' bands (Figures 
1 and 2). Our high-spatial resolution images achieved with AO reveal its complex morphology. 
The nebula has three bright clump-like structures at the following distances and position angles 
(P.A.) from LkHa 234: 2"78 ± 0"3 and 276° ± 4°, 2"78 ± 0"25 and 248° ± 1°, and 4"22 ± 
0"4 and 231° ± 1°. The nebula has the shape of a monopolar cavity, suggestive of the outflow 
activity frequently seen around Class I stars. A protostar might therefore be the source of 
the outflow. Infrared emission is often attributed to dust grains heated by a nearby source, 
which in this case may be an embedded protostar. However, the observed nebula extends more 
than 10". In addition, near-infrared polarization observations of the nebula showed that it had 
a typical concentric pattern explained by scattered light (Weintraub et al. 1994), indicating 
that it is a reflection nebula. The center of the polarization pattern does not coincide with 
LkHa 234, which also suggests that the nebula is not illuminated by LkHa 234, but by another 
object that is probably the source of the outflow. 

LkHa 234 is surrounded by another tenuous nebula that appears green around it in 



5 



Figure 3. This nebula extends approximately 5" (~6000 AU) toward the northeast from 
LkHa 234. No luminous stars can be seen except for LkHa 234 near this northeast neb- 
ula, suggesting that LkHa 234 probably illuminates it. Imaging polarimetry of this northeast 
nebula using the Lick AO system indicates its polarization position angle is consistent with 
its being scattered light from LkHa 234 (Perrin et al. 2006). The J — H and H — K colors of 
this northeast nebula measured at an offset (a, 5) = (— 0"43, +3"4) from LkHa 234 with an 
aperture radius of 20 pixels are 0.88 and 1.06 magnitudes, respectively. Meanwhile, the colors 
of LkHa 234 are 1.33 magnitudes and 1.12 magnitudes. Because the color of this nebula is more 
blue than that of LkHa 234, the scatterers of the starlight are not large grains (0.25-0.8 /jm) 
but probably smaller interstellar dust grains (< 0.25 /xm; Pendleton et al. 1990). 

Figure 4 shows the 11.8 /jm (SiC) and 17.65 jum images. These images show three dis- 
crete objects (LkHa 234, NW1, and NW2). NW1 and NW2 are located at 2"52 and 3'/69, 
respectively, from LkHa 234. These two objects are not seen in the J, if, and .fT-band im- 
ages; however, they are bright in the mid-infrared bands, suggesting that NW1 and NW2 are 
embedded young stars. In contrast, Objects F and G seen in the near-infrared images are 
not detected in the mid-infrared images, implying that they are relatively mature stars with 
less infrared excesses. Moreover, the 11.8 /mi image shows a diffuse envelope-like nebula east 
of LkHa 234 that curves northward of it. This is considered to be the same nebula as the 
tenuous northeast nebula discussed above. The region of diffuse 11.8 /im emission in Figure 
4 corresponds to the darker region just inside the northeast nebula's half-round structure in 
Figure 3. The curved mid-infrared nebula is probably a denser part of this northeast nebula. 
The mid-infrared nebular area just has a higher extinction; thus, it is opaque at less than 2 /mi, 
but it emits at longer wavelengths. 

4-2. YSO Candidates 

Eight YSO candidates (Objects B, C, D, E, F, G, NW1, and NW2) were detected at 
2" -11" from LkHa 234 in the near- and mid-infrared images. Table 2 lists the positions and 
photometric magnitudes of these objects. Table 3 lists objects observed by other authors that 
correspond to these objects. Four of them (Object D, E, F, and G) were detected in the K 
band for the first time. In particular, Objects F and G were discovered by our observations. 
Objects B and C correspond to IRS 4 and IRS 2, respectively named by Weintraub et al. (1994); 
NW1 corresponds to NW detected by Polomski et al. (2002) and IRS 6 detected by Cabrit et 
al. (1997). NW2 was detected at 11.8 and 17.65 /xm as well as in the V and M' bands. In total, 
nine stars of intermediate- and low-masses (LkHa 234 and Objects B, C, D, E, F, G, NW1, 
and NW2) were detected within 11" (~14000 AU) of LkHa 234, possibly forming a cluster. 
Below, we discuss the characteristics of the detected YSO candidates based on their positions in 
color-color diagrams, and we compare the present observational results with previous studies. 
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4-2.1. Color- color Diagrams 

Figure 5 shows JHK and JHKL' color-color diagrams. The infrared colors of five of the 
nine objects, LkHa 234, Objects B, C, D, and F, are plotted here. The three dotted lines denote 
the reddening law (Koornneef et al. 1983, Bessel & Brett 1988, Martin & Whittet 1990). If the 
color of an object is more red than that expected by interstellar extinction (the dotted lines), 
the object must have infrared excess in addition to its photospheric emission. Class I sources, 
which are in the stages of an accreting protostar (e.g., Lada 1987), are usually located to the 
right of the rightmost dotted line. Class II sources, which are objects like classical T Tauri 
stars with optically thick disks (e.g., Lada 1987), exist in the region between the rightmost and 
middle dotted lines. V band observations are useful for characterizing the infrared excesses of 
YSOs because deeply embedded objects can be more easily classified by their K — L' colors than 
their H — K colors (Lada et al. 2000). The JHKL 1 color-color diagram [Figure 5 (b)] indicates 
that LkHa 234 is more red than other objects (Objects B, D, and F), and it is located in the 
region where many Class I stars are located. For Object D, the lower limit of the J — H color 
is shown in Figure 5. The typical visual extinction toward the LkHa 234 region is Ay = 3.4 
(Hillenbrand et al. 1992). If Object D is a main-sequence star, its J — H color should be around 

3- 4 magnitudes and its visual extinction should be ~20 magnitudes based on the JHKL' color- 
color diagram. Because Object D is near the wall of the cavity of the reflection nebula, its high 
extinction may be because it is located behind the nebula. Therefore, Object D may be a YSO, 
but we cannot exclude the possibility that it is a main-sequence star. From the color-color 
diagrams, Objects B, C, and F are either main-sequence stars or Class III stars, which are 
objects like weak-lined T Tauri stars (e.g., Lada 1987). However, the fluxes of Objects B and C 
are significantly different than those from previous observations. We will discuss these objects 
in the next section. 

The stars we detect appear to be in distinct evolutionary phases from Class I to Class III 
or main-sequence, although they are located within the small region around LkHa 234. A 
similar characteristic has been observed in other star-forming regions. For example, Padgett 
et al. (2004) found that a region of 0.7 pc 2 in L1228 contains a total of five sources in various 
evolutionary stages. Note that the formation of NW1 and NW2 has occurred near LkHa 234, 
which is a relatively massive B5-7 type star, although the real separation between them remains 
unknown. 

4- 2.2. Variable Sources 

Objects B and C became brighter than the previous observations by Weintraub et al. 
(1994) in the H and K bands. The magnitudes of Object B decreased by 0.67 ± 0.24 magnitudes 
in the H band and 0.51 ± 0.16 magnitudes in the K band; those of Object C decreased by 
1.83 ± 0.23 magnitudes in the H band and 0.98 ± 0.16 magnitudes in the K band. Therefore, 
Object B and C are considered to be variable like typical YSOs. 

We now discuss the possible causes of this magnitude variation. First, we consider 
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the temporal variation of dust extinction. If the amount of the circumstellar and interstellar 
dust along the line of sight decreased after 1994, the object would become brighter. Adopting 
A H = 0.175 and A K = 0.112 when A v = 1 (Rieke & Lebofsky 1985), the increment of the H 
magnitudes Amu = +0.67 for Object B suggests that Ay changed by AAy = +3.82, correspond- 
ing to an increase in the iT-band magnitude of Amx = + 0.43, which is in good agreement with 
the observations. In contrast, Amu = +1.83 for Object C and the same calculation suggests 
AAy = + 10.46 corresponding to Amx ~ +1.17, which is not consistent with the observa- 
tions. Because the color-color diagrams indicate that Object B can be classified into Class II 
or Class III, there may be a circumstellar disk around Object B. Therefore, the variability of 
Object B may result from its obscuration by the inhomogeneous dust distribution in its disk. 
For Object C, an infrared excess has not been detected. However, the JHK color-color diagram 
alone cannot exclude the possibility of the existence of a disk (Lada & Adams 1992). 

Next, we consider the possibility that the variation occurs by the characteristics of the 
stellar surface. Cool spots are frequently observed for weak-line T Tauri stars, dwarfs, and 
the Sun. The temperatures of such cool spots are lower than the effective temperatures of 
their stars by ~1000 K (Stahler & Palla 2004), and these magnitude variations are caused by 
stellar rotation. To estimate the amplitude of magnitude variation for the case when a cool 
spot exists on the stellar surface, we used an equation from Stahler & Palla (2004): Am\ = 
—2.5 log[l — fs[l — B\(T spot )/ B\(T star )]]. The area ratio of the cool spot can be assumed to be 
20-40 percent (f B =0.2-0.4) of the stellar surface (Stahler & Palla 2004, Hatzes 1995, Vrba et 
al. 1988). In this case, the variation of magnitudes Am does not strongly depend on We 
derive the effective temperatures of Objects B and C from the J magnitudes by using the model 
of Baraffe et al. (1998; 2001). The age was assumed to be 10 6 yr. Although the estimated 
age of LkHa 234 is ~10 5 yr (Hillenbrand et al. 1992), Objects B and C may not be so young 
based on the color-color diagrams. If their age is 10 6 yr, their spectral types are F2-G2 (Baraffe 
et al. 1998; 2001). When these spectral types are assumed, the amplitudes of variation are 
~0.25 magnitudes in both the H and K bands. This variation is too small to account for the 
observed magnitude variations in Objects B and C; therefore, the cool spots are probably not 
responsible for these variations. 

Finally, we consider the possibility that hot spots cause magnitude variability. Hot spots 
are typically seen in classical T Tauri stars, and they are produced by the mass accretion from 
a circumstellar disk onto the stellar surface. The typical temperature of hot spots is ~7000 K 
if the stars are T Tauri stars, and hot spots cover on the order of 1 percent of the stellar 
surface area (Stahler & Palla 2004). We calculated the variation by using an equation from 
Stahler & Palla (2004): Am A ~ - 2.5 (log e)[(hc/ Xk B )(l/T star -l/T spot )\. For the case that the 
spectral types are F2-G2 and the temperatures of the hot spots are ~7000 K, the magnitude 
variations in the H band are Atuh = 1-5 ± 0.15 and in the K band are Am^ = 1.0 ± 0.1. This 
estimate demonstrates that hot spots can cause larger magnitude variations than cool spots. 
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Therefore, the variation seen in Object C may be caused by a hot spot on the stellar surface 
and not by circumstellar dust extinction because the hot spot model better reproduces the large 
amplitude of magnitude variation. Observations of periodicity will be required to discriminate 
between these two possibilities. In summary, the variation of Object B is likely caused by dust 
extinction, whereas that of Object C is likely caused by either dust extinction or hot spots. 

4-3. Jet-like Structures 

As shown in Figure 6, two jet-like structures are seen in the J and if-band images. These 
jet-like structures are located west of LkHa 234 with PAs of 253° ± 2° and 215° ± 2°. Because 
the [Fe II] emission line may contaminate the fluxes in the J and H bands, these structures 
may represent real jets such as those observed in L 1551 IRS 5 (Itoh et al. 2000, Pyo et al. 
2002). These two jet-like structures point neither to LkHa 234 nor NW1 but to Object G. 
The position of the northern jet-like structure is consistent with that of a [S II] jet detected by 
Ray et al. (1990), although the [S II] jet is more diffuse in their images perhaps partly due to 
lower spatial resolution caused by the seeing of l'.'b. [S II] and [Fe II] emissions usually have a 
similar molophology, and the emission is more compact in [Fe II] than in [S II] (e.g., Lopez et 
al. 2010). The PA of the [S II] jet was 252°, which is similar to that of our jet-like structure (PA 
= 253° ± 2°), but it is not clear whether this [S II] jet originates from Object G or LkHa 234 
because of the lower resolution. A H2 jet at PA = 226° was observed under the 0"5 seeing by 
Cabrit et al. (1997), and the position is roughly overlapped with that of the diffuse [S II] jet 
described above. However, the PA is different, and the H2 jet does not point to LkHct 234, 
Object G, or NW1, but to NW2. It is difficult to infer the driving source from the PA alone. 
Note that H2 and [Fe II] jets show different spatial distributions (e.g., Lopez et al. 2010), and 
jets sometimes bend and do not direct to the parent stars (e.g., Itoh et al. 2000). High-spatial 
resolution spectroscopy will help identify the driving sources of the jets (e.g., Pyo et al. 2002). 

4-4- Two Candidates for the Star Illuminating the Reflection Nebula 

A remarkable reflection nebula is seen west of LkHa 234. We consider two candidates 
besides LkHa 234 for the star illuminating the nebula: Object G and NW1. 

The position of Object G implies that it is the source of an outflow. First, it almost 
corresponds to the central position of the 2 /im polarization determined by Weintraub et al. 
(1994) (Table 4). In addition, the jet-like structures detected in the J and H bands point toward 
Object G. However, we are unable to exclude the possibility that Object G is a background 
star because it is difficult to estimate its luminosity using only its H- and i^-band magnitudes. 
Because the magnitudes in these bands may include the infrared excess, it is difficult to derive 
the amount of extinction, and thus its spectral type and distance. Therefore, it is also unclear if 
Object G is luminous enough to illuminate the reflection nebula in this region. Note, however, 
because Object G may be located near the mid-plane in the disk of NW1 or Object G itself, 
the extinction toward Object G may be quite high if it resides behind the disk mid-plane. If 
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Object G is associated with this region and affected only by the typical extinction toward this 
region (Ay =3.4, Hillenbrand et al. 1992), its mass is approximately that of a brown dwarf 
(with a spectral type of later than M6) based on evolutionary tracks (Baraffe et al. 1998; 2001, 
Siess et al. 2000). In this case, Object G might be the lowest- mass object found in this region. 
Previous studies have confirmed the existence of jets (~100 AU) and small-scale (~1000 AU) 
outflows driven by brown dwarfs and very low-mass stars (Phan-Bao et al. 2008, Lee et al. 
2009, Whelan et al. 2009). The jet-like structures and monopolar nebula we detected have 
larger sizes (> 1000 AU); therefore, if Object G is a brown dwarf, it might have difficulty 
driving the observed jets and outflow. In contrast, one might think that a higher extinction 
toward Object G is plausible as a high extinction of Ay = 15 is derived for Object F from 
the color-color diagrams (Figure 5). Assuming the same Ay as Object F, the spectral type of 
Object G is estimated to be M5 or K3 depending on the evolutionary tracks (Baraffe et al. 
1998; 2001, Siess et al. 2000). Therefore, we cannot exclude the possibility that Object G ejects 
these jets and outflow. 

NW1 is bright in the L', M', and mid-infrared bands (Figures 2 and 4), consistent with 
the conclusion by Fuente et al. (2001) that NW1 is an embedded young star. Figure 7 shows the 
spectral energy distribution (SED) of NW1; it was constructed using our near- and mid-infrared 
photometry and the millimeter flux densities at 1.3 and 2.6 mm obtained with the IRAM 30 m 
telescope (Fuente et al. 2001). Far- infrared data (60, 120, and 200 /mi) taken with the Infrared 
Space Observatory's (ISO) photometer (ISOPHOT) were also included in the SED (Abraham 
et al. 2000), but they are considered to be upper limits because ISO's large beam could not 
resolve NW 1 and NW 2. The dotted curve shows a gray-body sphere with a temperature of 
52 K and an outer radius of 667 AU. We adopted (3 = 1 and Tioo = 1 for r„= r^o where t wo 
is r at A = 100 /mi (O'linger et al. 2006, Tachihara et al. 2007). This curve is fitted to the 
millimeter data points, which represent the cold dust component (Fuente et al. 2001). The 
dashed curve shows the warm dust component with 100 K and 147 AU for the temperature 
and outer radius, respectively. Finally, the solid curve shows a sum of the cold and warm dust 
components. 

SED of NW1 that includes our new observations also indicates that NW1 is an embedded 
YSO. The mid-infrared and millimeter data can be fitted by combining the cold and warm dust 
components. The radius of the sphere derived from the fitting (667 AU) is consistent with 
the millimeter interferometric observations when the area of the continuum emission from cold 
dust grains is about 1000 AU in radius. In contrast, our near-infrared (V and M' bands) and 
mid-infrared (11.8 and 17.65 /mi) images show that the radius (FWHM) of NW1 is ~260 AU 
and ~140 AU, respectively. NW1 is extended in the V and M' bands, but scattered light 
may contribute to this; its MIR size better represents the thermal emission from warm dust. 
Therefore, the radius of the warm dust component inferred by the SED fitting (147 AU) is 
consistent with that infferred from our observed mid-infrared size (~140 AU). From the dual- 
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graybody fitting (solid curve), NWl's luminosity is estimated to be ~700 L & . The spectral 
type corresponding to this luminosity is B6-7 assuming it is a main-sequence star, which is 
similar to the B5-8 spectral type (~500 L Q ) suggested by Fuente et al. (2001). 

The following objects were also observed around LkHa 234: a redshifted CO outflow, a 
blueshifted [S II] jet, four radio continuum sources (VLA 1, VLA 2, VLA 3A, and VLA 3B), 
and associated H2O maser spots (Trinidad et al. 2004, Fuente et al. 2001, Umemoto et al. 2002; 
Table 4). H2O masers are usually found around accretion disks and outflows. As shown in 
Tables 2 and 4, the position of VLA 2 agrees well with that of NW2. NW2 is more faint in 
the near- and mid-infrared, implying that it is more embedded than NW1, or that it simply 
has a lower luminosity regardless of its evolutionary phase. Trinidad et al. (2004) favored 
the interpretation that VLA 2 is driving a large-scale CO and [S II] outflow based on the 
distribution of maser spots around it; however, the [S II] jet does not exactly point toward 
NW2 as seen in Figure 6. In contrast, NW1 is another candidate for a driving source of the CO 
outflow, as Fuente et al. (2001) originally suggested based on their interferometric CO J = 1 — >0 
images. VLA 3A and VLA 3B are located near the position of NW1. They are thought to 
form a close (~220 AU) binary system, and they seem to drive a thermal radio jet (Trinidad 
et al. 2004). NW1 appears as one source in our images, and its position is more consistent 
with that for VLA 3B. Moreover, NW1 is also located near the central position of the 2 /xm 
polarization (Weintraub et al. 1994; Table 4), and the morphology of the reflection nebula in 
the near-infrared strongly suggests that NW1 is associated with the monopolar nebula (Figures 
1 and 2). Consequently, it is possible that NW1 is a YSO of spectral type B6-7 that actively 
causes the large-scale CO outflow. If so, it could illuminate the reflection nebula in this region. 

In any case, LkHa 234 is not considered to drive the large-scale CO outflow and various 
jets (the [S II] and H 2 jets and the possible jets detected at J and H) observed in this region. 
We suggest that the source illuminating the reflection nebula is NW1, but we cannot completely 
exclude the possibility that it is Object G. 

4-5. Companion Candidate for LkHa 234 

Object F was identified for the first time. This object is located 1'.'9 (~2400 AU) from 
LkHa 234, and it is likely a star without infrared excess (see subsection 4.2). Assuming that 
the extinction Ay is 15 magnitudes (from the color-color diagrams of Figure 5) and the age is 
10 6 yr, we attempted to estimate the mass of Object F from its J-band magnitudes using model 
isochrones (Palla & Stahler 1999, Siess et al. 2000). We assumed an age of 10 6 yr for Object F 
since the age of LkHa 234 (10 5 yr) seems too young for the source without excess near-infrared 
emission. The estimate suggests that the spectral type of Object F is A2-A5 (2.0-2.5 M & ) if 
Palla & Stahler's (1999) method is adopted, whereas it is ~K4 if Siess et al's (2000) method 
is adopted. Therefore, Object F may be a low- or intermediate-mass star associated with this 
star-forming region, and a companion candidate for LkHa 234 separated from it by 2400 AU. 
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Herbig Ae/Be systems are frequently (30-75 percent) found in binary systems with separations 
between 30 and 4000 AU (Wheelwright et al. 2010, Carmona et al. 2007, Baines et al. 2006, 
Leinert et al. 1997). The mass ratio of LkHa 234 to Object F is 0.08-0.29, consistent with the 
results of previous observations of Herbig Ae/Be stars. If Object F is a companion, judging 
from the separation, it might have formed via disc fragmentation, which can produce a binary 
star with a higher mass ratio than that produced via core fragmentation (Wheelwright et al. 
2010, Kouwenhoven et al. 2009). 

The offset distance of Object G from LkHa 234, which is 2''3 (~2900 AU), is similar 
to that of Object F. If Object G is a companion of LkHa 234 and its mass is the same as 
that of a brown dwarf (Section 5.1), the mass ratio is below 0.01. This mass ratio is small 
compared with previous observations. In contrast, if the spectral type of Object G is ~M5 
or K3 (Section 5.1), the mass ratio is 0.03-0.1, which is slightly smaller than the mass ratio 
of Object F. There is little understanding about the complete census on the mass ratio and 
the separation because higher contrast is required to investigate fainter and closer companions 
around the Herbig Ae/Be stars (Wheelwright et al. 2010, Baines et al. 2006). In the case of 
Object G, it would be hard to determine if it is physically bound to LkHa 234, and to estimate 
the correct mass; however, our result demonstrates that AO imaging is useful for increasing the 
lower-mass sample of companion stars to Herbig Ae/Be stars to constrain the binary formation 
mechanism. 

5. Summary 

We observed the Herbig Ae/Be star LkHa 234 at near-infrared wavelengths (J, H, 
K, L' and M' bands) using CIAO on the Subaru Telescope. We obtained a high-angular 
resolution (0"2) with AO. The images reveal the complex morphology of the reflection nebula 
at approximately 3" from LkHa 234 that extends more than 10" westward. In addition, eight 
(Objects B, C, D, E, F, G, NW1, and NW2) YSO candidates were detected 2"-ll" from 
LkHo; 234, and four (Objects E, F, G, and NW2) of them were identified for the first time in 
the near-infrared. Object D is a Class II source or a main-sequence star, based on the JHK and 
JHKV color-color diagrams. Objects B and C are likely variable stars based on a comparison 
with previous photometry of them, and they could be YSOs. The variation of Object B is 
explained by variable extinction by circumstellar dust, and that of Object C may be caused by 
dust extinction or hot spots. LkHa 234, NW1, and NW2 were also detected in the 11.8 /xm 
(SiC) and 17.65 /im images obtained with the LWS on the Keck Telescope. NW1 and NW2 
are very bright in the mid-infrared bands, but they could not be seen in the J, H, and i^-band 
images. Therefore, they are likely embedded young stars. We found that these stars are in 
various evolutionary phases from Class I to Class III or main-sequence, and they are located 
within a small (0.14 pc 2 ) region around LkHa 234. Furthermore, we suggest that the source 
illuminating the reflection nebula is NW1 and not LkHa 234 although Object G may be another 

12 



candidate. NW1 is located near the central position of the 2 /im polarization and it is near the 
origin of the monopolar nebula. There is a possibility that NW1 is a YSO whose spectral type 
is B6-7 based on its SED and that it illuminates the reflection nebula in this region. Object F 
may be a companion star of LkHct 234 with a separation of 1'.'9 (~2400 AU) and a low-mass 
star with a spectral type of A2-A5 or ~K4 depending on the model isochrones. The mass ratio 
of LkHa 234 to Object F is then 0.08-0.29, consistent with previous observations of Herbig 
Ae/Be binaries. Moreover, Object G has a similar separation from LkH« 234 (~2900 AU), and 
if it is a companion, the mass contrast may be even larger. 

We thank D. Kontopoulos for the careful reading of our manuscript. We are also grateful 
to T. Kohyama for the kind support. 
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Fig. 2. (a) L'-band image (left panel) and (b) M'-band image (right panel) of the LkHce 234 region. These 
images are expressed in logarithmic scale. North is up and east is to the left. Five objects (Objects D, E, 
F, NW1, and NW2) are detected around LkHa 234. NW1 is bright in these bands. The display level of 
(a) is 1 a to 8.5 a, while that of (b) is 1 a to 5.7 a. 







Table 1. 


Summary of 


near- 


infrared observations 






Date (yy-mm-dd) 


Band 


T eX p [s] 


Coadd 


It ot 


[s] Mask size ["] 


Standard star 


Seeing ["] 


06-07-15 


J 


10 


1 


60 


0.8 


FS 150 


0.8 




H 


6 


1 


300 


0.8 


FS 150 


0.8 




K 


6 


1 


300 


0.8 


FS 150 


0.8 


08-05-25 


J 


10 


1 


300 


0.8 


FS 150 


0.5 


08-07-10 


V 


0.18 


10 


252 




HD129655 


1.3 




M' 


0.18 


10 


468 




HD129655 


1.3 


04-07-26 


11.8 


0.01 


28500 


285 




HD197989 




04-08-28 


17.65 


0.01 


28500 


285 




HD197989 





Columns 1 and 2 list the dates (year-month-date) and the observational bands. Columns 3 and 4 represent 
the exposure time and the number of coadds for each frame. T to t in Column 5 gives the total exposure time. 
Column 6 gives the diameter of the occulting mask; we utilized a circular Lyot stop at the pupil plane to 
block out 20 percent of the pupil diameter. Column 7 lists the photometric standard stars (Leggett et al. 
2003, Leggett et al. 2006, and Cohen et al. 1999). Column 8 gives the natural seeing in optical. 
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Fig. 3. JHK image displayed in logarithmic scale to show the faint northeast nebula (indicated by the 
ellipse). The faint nebula does not appear in Figure 1 because of the different display scale. The contour 
level is 0.7 a to 2.7 a in the ff-band image. Ghost images, like the one in northeast of Object C, are 
masked in this figure. Note the half-round structure of this nebula at approximately 5" (^6000 AU) 
northeast of LkHa 234. North is up and east is to the left. 
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Fig. 4. (a) 11.8 /im (SiC; left panel) and (b) 17.65 /im (right panel) Keck LWS images around LkHa 234 
displayed in logarithmic scale. North is up and east is to the left. The field of view is 9" x 9" for both. 
Contours are shown logarithmically spaced from 0.001 to 0.33 of the peak intensity for the 11.8 /xm image 
and from 0.01 to 0.33 of the peak intensity for the 17.65 /on image. Note that the embedded sources of 
NW1 and NW2 are bright at the mid-infrared wavelengths. In particular, NW1 is comparable in flux to 
LkHa 234 at 17.65 /im. 



Table 2. Positions and magnitudes of detected objects. 



Object 


Dis. [" 


]p.a. [°; 


J H K 


V M' 


SiC [Jy] 17.65 [Jy] 


LkHa 234 


: 0.0 


0.0 


9.528 (0.02)8.201 (0.02)7.081 


(0 


.02)5.23 (0.05) 4.83 (0.18) 3 


.53 (0.15)4.95 (0.50) 


B 


8.28 


11.16 


15.05 (0.06)13.43 (0.04)12.39 


(0 


.04)11.45 (0.06)11.19 (0.36) 




C 


11.12 


265.3 


13.15 (0.06)12.25 (0.03)11.88 


(0 


.04) 




D 


10.70 


205.4 


>18.5 15.50(0.10)13.17 


(o 


.06)11.40 (0.05)11.10 (0.56) 




E 


6.37 


204.4 


>18.5 >19.0 15.04 


(o 


.30)10.39 (0.05)9.40 (0.22) 




F 


1.89 


96.9 


15.45 (0.25)13.29 (0.23)11.99 


(0 


.22)11.6 (0.06) >11.37 




G 


2.28 


275.9 


>18.5 18.26 (0.40)16.46 


(o 


.20) 




NW1 


2.52 


295.7 






9.50 (0.06) 8.32 (0.2) 


.38 (0.02)3.64 (0.40) 


NW2 


3.69 


305.5 






10.58 (0.06)9.56 (0.4) 


.03 (0.01)0.95 (0.30) 



The distances and PAs are from LkHa 234. The uncertainty in the offset distance and in the position angle at 
the near-infrared wavelengths are 0."05 and l.°7, respectively. Those at the mid-infrared wavelengths are 0."15 and 
2.°6, respectively. The figures in parentheses are the uncertainties (±) in the photometric values. The detection 
limits in the J, H, K, L' , and M' bands are 18.5 ± 0.2 magnitudes, 19.0 ± 0.3 magnitudes, 17. 7± 0.2 magnitudes, 
12.1 ± 0.2 magnitudes, and 10.0 ± 0.2 magnitudes, respectively. The magnitudes of LkHa 234 in the JHK bands 
are from the 2MASS All-Sky Catalog. 
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Fig. 5. (a) JHK color-color diagram (left panel), (b) JHKL' color-color diagram (right panel). Five 
objects (Objects B, C, D, F and LkHa 234) are plotted. The JHK colors are transformed to the CIT 
photometric system (Cutri et al. 2003, Itoh et al. 2001). The bold curve shows colors of unrcddcncd 
main-sequence (Lada & Adams 1992 and Koornneef et al. 1983) and T Tauri (Meyer et al. 1997) stars. 
The intrinsic colors of the main sequence were determined from the values for V — K as a function of the 
spectral type (Koornneef et al. 1983), and spectral types from 06 to M8 were included. The three dotted 
lines are parallel to the reddening vector. In the JHK color-color diagrams, the gradient of the reddening 
vector is 1.7 (Koornneef et al. 1983); in the JHKL 1 color-color diagram it is 2.09 assuming the reddening 
law of Bessel & Brett (1988) and Martin & Whittet (1990). Crosses represent the amounts of extinction 
(Ay), which are 0, 5, 10, 15, 20, 30 and 40 magnitudes. 



Table 3. Corresponding objects in previous studies. 



Near/Mid IR objects (in this paper) 


Corresponding object 


Object B 


IRS 2 (e.g., Weintraub et al. 1994) 


Object B 


IRS 4 (e.g., Weintraub et al. 1994) 


Object E 


C3 (Tofani et al. 1995) 


LkHa 234 NW1 


PS1 (Weintraub et al. 1994) 




10 fxm companion or IRS 6 




(Cabrit et al. 1997, Fuente et al. 2001) 




VLA 3A+B (Trinidad et al. 2004) 




LkHa 234 NW (Polomski et al. 2002) 


LkHa 234 NW2 


VLA 2 (Tofanil995) 



20 



LkHa 234 





Fig. 6. (a) J-band and (b) ff-band images around LkHa 234. North is up and east is to the left. The 
white dashed arrows indicate two jet-like structures located west of LkHa 234. The PAs of these structures 
arc 251°-255° and 214°-217°. The contour levels of these J and H images are 64 a to 128 a and 51.5 a 
to 110 a, respectively. Blue and red arrows represent [S II] and H2 jets, respectively (Ray et al. 1990, 
Cabrit et al. 1997). Crosses indicate Object G, NW1 and NW2. Filled circles indicate VLA 2, VLA 3A 
and VLA 3B. A filled triangle indicates the mid-infrared object IRS 6 discovered by Cabrit et al. (1997). 
The two jet-like structures point to neither LkHa 234 nor NW1 but to Object G. The position of the 
[S II] jet is roughly overlapped with those of the jet-like structures and H2 jet, although the [S II] jet is 
very diffuse. The H 2 jet does not point to LkHa 234, Object G, or NW1, but to NW2. 



Table 4. Positions of outflow source candidates from other observations. 



Object 


Distance ["] 


P.A. [°] 


Author 


PS 1 


3 


290-340 


Weintraub ct al. 1994 


PS 1 


2.58 




Weintraub et al. 1996 


10 [ira companion 


2.7 


280-286 


Cabrit et al. 1997 


VLA 1 


5.91 


312 


Tofani et al. 1995 


VLA 2 


3.50 


300 


Tofani et al. 1995 


VLA 3A 


2.23 


292 


Trinidad ct al. 2004 


VLA 3B 


2.40 


295 


Trinidad et al. 2004 



Offset distances and offset position angles are from LkHa 234. 
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Fig. 7. Spectral energy distribution of NW1. The mid-infrared (11.8 /an and 17.65 /im) data are from 
the present study. The far-infrared (60, 120 and 200 /im) and millimeter (1.3 and 2.6 mm) data are from 
Abraham et al. (2000) and Fuente et al. (2001), respectively. The far-infrared observations are only upper 
limits because NW1 is unresolved. The dotted curve represents a gray-body sphere with temperature of 
52 K and an outer radius of 667 AU; this curve is fitted to the millimeter data observations (cold dust 
components). The dashed curve represents warm dust components (at 100 K and 147 AU). The solid 
curve shows the result of adding the cold and warm dust components. These fittings were adopted using 
/3 = 1 and t wq = 1 (O'linger et al. 2006 and Tachihara et al. 2007). 
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